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ABSTRACT: The evaporation coefficient and equilibrium vapor pressure of silicon 
monoxide over a mixture of silicon and vitreous silica have been studied over the 
temperature range (1433 to 1608 ) K. The evaporation coefficient for this temperature 
range was ( 0.007 ± 0 . 002 ) and is approximately an order of magnitude lower than the 
evaporation coefficient over amorphous silicon monoxide powder and in general 
agreement with previous measurements of this quantity. The enthalpy of reaction at 
298.15 K for this reaction was calculated via second and third law analyses as (355 ± 25 ) 
kJ-moT 1 and ( 363.6 ± 4 . 1 ) kJ-moU 1 , respectively. In comparison with previous work with 
the evaporation of amorphous silicon monoxide powder as well as other experimental 
measurements of the vapor pressure of silicon monoxide gas over mixtures of silicon and 
silica, these systems all tend to give similar equilibrium vapor pressures when the 
evaporation coefficient is correctly taken into account. This provides further evidence that 
amorphous silicon monoxide is an intimate mixture of small domains of silicon and silica 
and not strictly a true compound. 
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20 ■ INTRODUCTION 

21 The thermodynamics of the silicon— oxygen (Si— O) system are 

22 extremely important in understanding and improving the 

23 growth of silicon monoscrystals grown by the Czochralski 

24 (Cz) method for the semiconductor industry. In these systems 

25 the transport of oxygen in the silicon melt to the crystal growth 

26 interface plays an important role in the final properties of the 

27 resulting silicon wafers . 1 One of the challenges of growing 

28 larger crystals via the Czochralski process is the attack of the 

29 molten silicon on the silica glass crucible that occurs during the 

30 longer process times required for these larger crystals . 2 ' 3 These 

31 pits form as small particles break off from the crucible and may 

32 be transported to the crystal interface, forming dislocations. 

33 Another potential problem is that the silicon monoxide (SiO) 

34 gas liberated during the crystal growth may condense above the 

35 melt, forming deposits that could potentially fall back into the 

3 6 melt and also contaminate the growing crystal . 3 Data on the 

37 silicon— silica system and the vapor pressure of SiO over such a 

38 mixture can be useful in modeling and optimizing Cz— Si 

39 growth. 4 5 

40 Another area in which the properties of the Si— O system are 

41 important is the modeling of grain formation in stellar outflows. 

42 One of the most abundant species in the outflows of oxygen- 

43 rich, asymptoptic giant branch (AGB) stars is believed to be 

44 silicon monoxide. Despite its relatively high abundance, 

45 previous models of these outflows seem to indicate that silicon 

46 monoxide would not appreciably nucleate until approximately 

47 600 K, well below the observed formation of silicate grains 

48 above 1000 K in these oxygen-rich stars . 6 Therefore, other less 

49 abundant, but more refractory species were theorized to form a 


seed nucleus upon which silicon monoxide could later so 
condense. si 

In a previous work, the vapor pressure and evaporation 52 
coefficient of silicon monoxide over amorphous silicon 53 
monoxide powder was measured . 7 These data coupled with 54 
previous measurements of SiO vapor pressure showed that the 55 
actual vapor pressure of SiO at approximately 1000 K and 56 
below was much lower than the expression that had been 57 
previously used in modeling SiO nucleation and growth , 8 ' 9 58 
thereby greatly increasing the possibility that SiO is the initial 59 
condensate in oxygen-rich, stellar outflows . 10 so 

In this earlier work, both the evaporation coefficient and the 61 
vapor pressure of SiO( g j over SiO^) were measured. In this 62 
present study, these same quantities are measured for the 63 
following reaction: 64 

1 / 2Si (cr) + l / 2 Si0 2 ( vit ) -*■ SiO(g) (!) 

In both the case of this silicon and silica reaction as well as the 65 
sublimation of amorphous silicon monoxide, the product is 66 
silicon monoxide gas. Brewer and Edwards reviewed the 67 
available thermodynamic and spectroscopic data on the 68 
silicon— oxygen system and concluded that the reaction given 69 
by eq 1 yields virtually pure silicon monoxide gas and that this 70 
gas is a monomer . 11 Later, this thermodynamic assessment was 71 
experimentally proven by Porter et al . 12 Using mass 72 
spectrometry, these authors were able to verify that the 73 
predominant species in the vapor over eq 1 was silicon 74 
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75 monoxide gas, with Si 2 0 2 being the next abundant species at a 

76 concentration approximately 4 orders of magnitude smaller. 

77 For many years, there has been a controversy as to whether 

78 SiO^) is a compound or simply an intimate, stoichiometric 
7 ? mixture of silicon (Si) and silica (Si 0 2 ). 13 Recent works seem 
so to confirm that SiO( am ) is actually a mixture with very fine 

81 domains (~5 nm) of Si and Si 0 2 . 4,15 In studies of the vapor 

82 pressure of silicon monoxide gas, low evaporation coefficients 

83 have been noted for both systems. For example, the 

84 evaporation coefficient for SiO^) has been measured as 

85 approximately 0.05 while the evaporation coefficient for SiO( g ) 

86 over a mixture of Si and Si 0 2 is typically an order of magnitude 

87 smaller. Such low evaporation coefficients were measured by 

88 Rocabois et al. for both systems . 16 Furthermore, these authors 

89 note that, within the experimental uncertainty for their system, 

90 both systems yielded the same vapor pressure for silicon 

91 monoxide gas when the evaporation coefficients were taken 

92 into account. In a previous work, both the evaporation 

93 coefficient and the vapor pressures for silicon monoxide gas 

94 over SiO( am ) were in very good agreement with the results from 

95 Rocabois et al . 7 In this work, we report measured evaporation 

96 coefficients and equilibrium vapor pressures for silicon 

97 monoxide gas over a mixture of silicon and silica and compare 

98 these data with previous measurements. 

99 ■ EXPERIMENTAL APPARATUS AND PROCEDURE 

100 The vapor pressure over the silicon— silica mixture is measured 

101 by monitoring the mass loss from a Knudsen effusion cell using 

102 a thermogravimetric balance. The thermogravimetric balance 

103 used in this work, a Thermo-Cahn 2171 , is capable of reaching 

104 temperatures of 1973 K. The furnace tube of this balance is 
ios connected to mechanical and turbomolecular pumping system 

106 that maintains the pressure < 1 (T 3 Pa in the sample area. The 

107 resulting mass loss from the sample cell is continuously 
ios recorded with a computer connected to the balance. During an 
109 experimental run, the sample cell is suspended on the sample 
no side of the balance arm and centered within the furnace section 
in of the reaction tube. This section is heated using six resistive 

112 heating elements. The temperature within the furnace is 

113 measured with a type-B thermocouple sheathed in an alumina 

114 tube and positioned approximately 0.5 cm below the bottom of 
ns the Knudsen cell. 

116 The Knudsen cell used in this work is constructed from the 

117 closed-end sections of two different-sized, 99.8 % alumina tubes 
ns as shown in Figure 1 . A smaller inner tube with a drilled orifice 

119 is cut and inverted in a close-fitting, larger tube. A window in 

120 the outer tube is cut, and the orifice of the inner tube is 

121 positioned in this this window. Although there is little clearance 

122 between the tubes, the tubes are also sealed at the top and the 

123 window areas with a zirconia-based adhesive (Resbond 904 , 

124 Cotronics Corp.). All parts in the hot furnace zone area are 

125 made of 902 alumina with the exception of this zirconia 

126 adhesive, and in this case the minimum amount necessary to 

127 seal the cell is used. 

128 Diamond-coated drill bits are used to produce the small 

129 effusion orifices in the alumina cells. Because of their small size, 

130 very high speed rotational rates are needed, and the bit must be 

131 fed at extremely slow feed rates. In this work a high speed air 

132 grinder (~ 8 S 000 rpm) was used, and it was fed using a 
•.33 computer controlled stage. This allowed the production of very 

134 clean orifices down to approximately 0.5 mm. 

13 5 Before any experiments were performed, an equilmolar 

136 mixture of silicon powder (100 mesh, Alfa Products) and 


Anic'c 



Figure 1 . Alumina efiusion cell constructed from two closed-end 
tubes. 


silicon dioxide (Mathey "Specpure” grade) was prepared and 137 
stored in a sealed container. Samples for all of the runs in this 138 
work were taken from this same container. For an experiment, 139 
alumina tubes were cut to construct a cell as shown in Figure 1 , 140 
and an orifice was then drilled in the smaller tube. To measure 141 
the orifice diameter, the cell was placed under a microscope and 142 
measured with the aid of a traveling stage of micrometer 143 
accuracy . 7,17 After filling the cell with the Si/Si 0 2 mixture, the 144 
cell was sealed with the zirconia cement and allowed to cure. 14s 
After curing, the cell would be placed on the sample side of the 146 
balance within the furnace tube and evacuated for 24 h at room 147 
temperature. The sample was then heated to various temper- 14s 
atures, and the mass of the sample cell was continuously 149 
monitored, taking temperature and mass measurements twice 150 
per second. 151 

In the same manner as was done for a previous study with 152 
silicon monoxide, after a run of 2 min, averages of the resulting 153 
large data set were taken for both the temperature and the mass 154 
data points. These data were then processed by a computer 155 
program to calculate the mass loss rates as a function of 156 
temperature. In this program, the temperature data were 157 
searched for regions where the temperature remained 158 
essentially constant. For these isothermal periods, the rates of 159 
mass loss were constructed from these 2 min averages. All of 160 
these values for an isothermal period were then averaged, and 161 
the standard deviation of this mean value was then used as an 1 62 
estimate of the uncertainty in this mass loss rate and used in the 163 
calculation of the vapor pressure. 164 

■ VAPOR PRESSURE CALCULATION i*s 

The measured vapor pressure, P ^ is related to the mass loss 166 
rate, m, via the Hertz— Knudsen equation: 167 

p _ m j2nRT 

(2) 

where B is the cross-sectional area of the effusion orifice, R is 168 
the gas constant, T is the temperature of the gas, and M w is the 169 
molecular weight of the effusate. The term, W B , is the Clausing 170 
correction factor for the cell orifice. If the cell wall thickness is 171 
knife-edge thin, this factor is 1 . 0 . If, as is typical, the wall in the 172 
vicinity of the orifice has some finite thickness, this short "pipe” 173 
can cause an impedance to the transmission of molecules from 174 
the cell. This factor, W Bt accounts for this back-reflection of its 
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176 some of the molecules and is based on the work of Clausing . 18 

177 For arbitrary geometries, this factor must be computed from 

178 the Clausing integral equation. Berman developed a series 

179 expansion approximation to this equation for capillaries, and 
iso this formula has been used to calculate the Clausing factors for 
isi the cell . 19 For an effusion orifice of radius, a, and wall thickness 

182 (“pipe" length), l, the Clausing coefficient, W, is given by 

W *Ql-Q .2 ( 3 ) 

183 where 

Q.! = 1 + (L 2 / 4 ) - (L/ 4 )(L 2 + 4) 1/2 ( 4 ) 


d 2 = [(8 - L 2 )(L 2 + 4) 1/2 + L 3 - 16] 2 
/[72 L(L 2 + 4) 1/2 

- 288 In[L + (L 4 + 4 ) 1/2 ] + 288 In 2 ] ( 5 ) 

284 and 


L = \/a (6) 

185 The accuracy of this equation has been verified by Monte Carlo 

186 simulations , 20 and this expression is reported to be better than 

187 0.1 % for 0 < L < 5 , a range common to Knudsen effusion 

188 studies . 21 

! 89 As noted in the Introduction, silicon monoxide is expected to 

190 have a low evaporation coefficient. Whitman and Motzfeldt 

191 developed the following equation to account for such low 

192 evaporation coefficients in Knudsen cells : 22>23 



193 In this equation the impedance of the flow of molecules to the 

194 effusion orifice height is included in the term given by W A , the 

195 Clausing factor for the cell. The term, P^, is the true, 

196 equilibrium vapor pressure, a, is the evaporation coefficient, 

197 and the term / is a factor related to the cell geometry and is 

198 given by 

1 Ak (8) 


199 where A is the cross-sectional area of the cell and k is the ratio 

200 of the effective evaporation surface area to the cell cross- 

201 sectional area. Therefore, the product Ak gives the actual 

202 evaporation surface area. In general, k is very difficult to 

203 accurately quantify and will depend upon a variety of factors 

204 including the packing of the sample, how finely divided the 

205 sample is, and so forth. In this current work, the value of k is 

206 assumed to be 1 . 

20 r In experiments with low evaporation coefficients, this value, 
208 a, can be estimated by rearranging eq 7 to 


p m 


= p eq - pj\ 



( 9 ) 


209 Therefore, if a series of effusion experiments are made with 

210 cells of different orifice sizes and the data were plotted on a 

211 graph as P m versus P„J and fit to a straight line, the resulting 

212 intercept should yield the equilibrium vapor pressure, A* while 

213 the value of the evaporation coefficient, a, may be computed 

214 from the slope. 


A few of additional notes should be made regarding eq 9 . 2:5 
First, for orthocylindrical cells, the value of W A is 0 . 5 , and the 216 
quantity ((l /Wj) — 2 ) = 0 . In this work, the value of W A is 217 
close to 0.5 so the sum of these two terms contribute very little 218 
to the sum in the parentheses of eq 9 , especially for cases where 219 
a is very small. Next, as mentioned earlier the value of k was 220 
assumed to be 1 . For values of k > 1 , the value of a will be 221 
smaller, and the evaporation coefficients reported in this work 222 
are, strictly speaking, ka products, and this should be noted 223 
when comparing evaporation coefficient data from other 224 
experimental works. 225 

■ ESTIMATION OF EXPERIMENTAL UNCERTAINTIES 226 

The estimated uncertainties in the current work are similar to 227 
those for a previous work with silicon monoxide, and the reader 228 
is referred to this work for more detail . 7 As noted in the 229 
previous section, values of the equilibrium vapor pressure and 230 
evaporation coefficient are taken from plots of P m versus at 231 
particular temperatures for different effusion orifices. Un- 232 
certainties in these quantities arise from the measured 233 
temperature, the measured mass loss rate, and cell geometry 234 
(cell wall thickness, orifice diameter). 23s 

Typical uncertainties for a type-B thermocouple in the 236 
experimental temperature range are ± 5 K, and these values are 237 
used in this work. As in the case of experimental runs with 238 
silicon monoxide, the accuracy of the thermocouple was 239 
checked against the melting point of a sample of pure copper 240 
with the thermocouple falling within this ± 5 K-band of the 241 
copper fusion temperature . 7 242 

Factors related to cell geometry include the orifice diameter, 243 
the depth of the effusion orifice, and the interior cell geometry. 244 
The effusion orifice was measured using a microscope and an 245 
accurate, moveable stage, and these measurements were 246 
estimated to be accurate to within ± 0.02 mm. The wall 247 
thickness at the effusion orifice was constant for all the cells 248 
used and was measured with a similar uncertainty as ( 1.27 ± 249 
0 . 02 ) mm. To calculate the Clausing factor for the effusion cell 250 
used in the Whitmann— Motzfeldt equation, the distance from 251 
the top of the evaporating material to the effusion orifice must 252 
be known. These interior measurements are not accurately 253 
known. Fortunately, the results are not sensitive to these values, 254 
especially in the case of very low evaporation coefficients. In 255 
this work as well as in previous experiments with silicon 256 
monoxide powder, the height from the top of the evaporating 257 
material to the effusion orifice was estimated as (10 ± 2 ) mm, 2S8 
while the evaporating surface area was just taken as the internal 259 
cross-sectional area of the inner tube with a diameter of ( 6.3 ± 260 
0 . 2 ) mm. 261 

For an isothermal period, the mean mass loss rates of the 2 262 
min averages described earlier are used as the mass loss rate in 263 
eq 2 , and the uncertainty in this quantity is estimated as the 264 
standard deviation of the measured value. Typically, the 265 
predominant source of uncertainty in the calculation of P m in 266 
eq 2 is the mass loss rate. Nevertheless, the uncertainties in all 267 
of these measured observables, (temperature, mass loss rate, 268 
and cell geometry), on P ^ are considered, and these estimates 269 
are added in quadrature as described in the earlier publication 270 
for silicon monoxide. 271 

■ RESULTS 272 

The rate of mass loss from several different effusion cells was 273 
monitored over time at several specific temperatures. After- 274 
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Table 1 . Experimental Data for SiO Vapor Pressure Measurements over a Mixture of Silicon and Silica 1 * 


Article 


r 

d 

B 


tn 

r 

m 


P m 

k 

mm 

cm 2 

Wb 

mg 

s 

rngmiri -1 

/ 

Pa 

1433 

1.08 

0.00925 

0.479 

15.34 

21000 

0.0439 ± 0.0050 

0.0140 

2.15 ± 0.27 

1433 

1.24 

0.01217 

0.510 

18.64 

210 C 0 

0.0532 ± 0,0034 

0,0196 

1.86 ± 0.15 

1433 

0.57 

0.00257 

0.339 

6.82 

32100 

0.013 ± 0.011 

0.0027 

3.2 ± 2.9 

1433 

1.72 

0.02324 

0,583 

21,00 

21300 

0.0592 ± 0.0057 

0.0428 

0.950 ± 0.096 

1443 

1.72 

0.02324 

0.583 

21.30 

17400 

0.07346 ± 0.00058 

0.0428 

1.183 ± 0.036 

1443 

0.57 

0.00257 

0.339 

6.30 

17100 

0.02209 ± 0.00073 

0,0027 

5.53 ± 0.99 

1443 

1.08 

0.00925 

0.479 

16.02 

17400 

0.05522 ± 0.00047 

0.0140 

2.72 ±0.16 

1443 

1.24 

0.01217 

0.510 

17.97 

17400 

0.06195 ± 0.00069 

0.0196 

2.18 ± 0.10 

1453 

0.57 

0.00257 

0.339 

6.33 

14700 

0.02585 ± 0.00078 

0.0027 

6.5 ± 1.2 

1453 

1.72 

0.02324 

0.583 

21.15 

14700 

0.08636 ± 0.00062 

0.0428 

1.395 ± 0.043 

1453 

1.08 

0.00925 

0.479 

15.66 

14700 

0.0639 ± 0.0011 

0.0140 

3,16 ± 0.19 

1453 

1.24 

0.01217 

0.510 

17.36 

14700 

0.07084 ± 0.00069 

0.0196 

2.50 ± 0.12 

1442 

1.72 

0.02324 

0.583 

21.31 

12600 

0.10152 ± 0.00061 

0.0428 

1.645 ± 0.050 

1442 

1.24 

0,01217 

0.510 

17.33 

12600 

0.08253 ± 0.00069 

0.0196 

2.92 ± 0.14 

1462 

1.08 

0.00925 

0.479 

15.71 

12600 

0.0749 ± 0.0020 

0.0140 

3.71 ± 0.23 

1462 

0.57 

0.00257 

0.339 

6.28 

12300 

0.0307 ± 0.0011 

0.0027 

7.7 ± 1.4 

1472 

1.08 

0.00925 

0.479 

15.64 

1050 C 

0.0893 ± 0.0035 

0.0140 

4.44 ± 0.31 

1472 

1.24 

0.01217 

0.510 

16.73 

10500 

0.09560 ± 0.00079 

0.0196 

3.39 ± 0.16 

1472 

1.72 

0.02324 

0.583 

20.91 

10500 

0.11949 ± 0,00069 

0.0428 

1.943 ± 0.059 

1472 

0.57 

0.00257 

0.339 

6.57 

10500 

0.0376 ± 0.0012 

0.0027 

9.5 ± 1.7 

1431 

0.57 

0.00257 

0.339 

6,77 

9000 

0.0452 ± 0.0079 

0.0027 

11.5 ± 2.8 

1431 

1.72 

0.02324 

0.583 

21.09 

9000 

0.14060 ± 0.00072 

0.0428 

2.293 ± 0.069 

1431 

1.08 

0.00925 

0.479 

15.58 

9000 

0.10391 ± 0.00093 

0.0140 

5.19 ± 0.30 

1431 

1.24 

0.01217 

0,510 

16.96 

9000 

0.11307 ± 0.00095 

0.0196 

4.02 ± 0.19 

1491 

1.08 

0.00925 

0.479 

15.53 

7500 

0.1242 ± 0.0012 

0,0140 

6.22 ± 0.36 

1491 

1.24 

0.01217 

0.510 

16.77 

7500 

0.1341 ± 0.0012 

0.0196 

4.79 ± 0.23 

1491 

1.72 

0.02324 

0.583 

21.40 

7800 

0.16467 ± 0.00091 

0.0428 

2.694 ± 0.081 

1491 

0.57 

0.00257 

0.339 

7.38 

7500 

0.0590 ± 0.0014 

0.0027 

15.0 ± 2.7 

1500 

1.08 

0.00925 

0.479 

15.54 

6300 

0.14790 ± 0.00079 

0.0140 

7.43 ± 0.42 

1500 

1.24 

0.01217 

0.510 

16.48 

6300 

0.1569 ± 0.0010 

0.0196 

5.62 ± 0.26 

1500 

1.72 

0.02324 

0.583 

21.12 

6600 

0.1921 ± 0.0013 

0.0428 

3.154 ± 0.096 

1500 

0.57 

0.00257 

0.339 

6.95 

6300 

0.0661 ± 0.0016 

0.0027 

16.9 ± 3.0 

1510 

1.72 

0.02324 

0.583 

20.04 

5400 

0.2225 ± 0.0016 

0.0428 

3.66 ± 0.11 

1510 

1.08 

0.00925 

0.479 

15.44 

5400 

0.1717 ± 0.0013 

0.0140 

8.65 ± 0.50 

1510 

1.24 

0.01217 

0.510 

16.56 

5400 

0.18399 ± 0.00093 

0.0196 

6.61 ± 0.31 

1510 

0.57 

0.00257 

0.339 

7.01 

5400 

0.07790 ± 0.00090 

0.0027 

19.9 ± 3.5 

1519 

1.24 

0.01217 

0.510 

16.10 

4500 

0.2146 ± 0.0011 

0.0196 

7.74 ± 0.36 

1519 

0.57 

0.00257 

0.339 

7.16 

4800 

0.0893 ± 0.0017 

0.0027 

22.9 ± 4.1 

1519 

1.08 

0.00925 

0.479 

14.96 

4500 

0.19950 ± 0,00098 

0.0140 

10.08 ± 0. S 8 

1520 

1.72 

0.02324 

0.583 

20.52 

4800 

0.2567 ± 0,0018 

0.0428 

4.24 ± 0.13 

1530 

1.72 

0.02324 

0.583 

19.63 

3900 

0.3018 ± 0.0017 

0.0428 

5.00 ± 0.15 

1530 

1.08 

0.00925 

0.479 

15.36 

3900 

0.2363 ± 0.0011 

0.0140 

11.98 ± 0.68 

1530 

1.24 

0.01217 

0.510 

16.53 

3900 

0.2543 ± 0.0014 

0.0196 

9.20 ± 0.43 

1530 

0.57 

0.00257 

0.339 

6.79 

3900 

0.1043 ± 0,0016 

0.0027 

26.9 ± 4.7 

1540 

0.57 

0.00257 

0.339 

6.35 

3300 

0.1151 ± 0.0020 

0.0027 

29.8 ± 5.3 

1540 

1.72 

0.02324 

0,583 

19.02 

3300 

0.3457 ± 0.0022 

0.0428 

5.75 ± 0.17 

1540 

1.08 

0,00925 

0.479 

14.99 

3300 

0.2725 ± 0.0019 

0.0140 

13.86 ± 0.79 

1540 

1.24 

0,01217 

0.510 

16.21 

3300 

0.2946 ± 0,0020 

0.0196 

10.69 ± 0.50 

1549 

0.57 

0.00257 

0.339 

6.32 

2700 

0.1398 ± 0.0046 

0.0027 

36.3 ± 6.5 

1549 

1.72 

0.02324 

0.583 

17.89 

2700 

0,3972 ± 0.0023 

0.0428 

6.63 ± 0.20 

1549 

1.08 

0.00925 

0.479 

14.08 

2700 

0,3127 ± 0.0021 

0.0140 

15.96 ± 0.91 

1549 

1.24 

0.01217 

0.510 

15.28 

2700 

0,3397 ± 0.0023 

0.0196 

12.37 ± 0.58 

1560 

0.57 

0.00257 

0.339 

4.64 

1800 

0.1540 ± 0.0032 

0.0027 

40,1 ± 7.1 

1560 

1.08 

0.00925 

0.479 

12.78 

2100 

0.3648 ± 0,0026 

0.0140 

18,7 ± 1.1 

1560 

1.24 

0.01217 

0.510 

13.95 

2100 

0.3977 ± 0.0036 

0.0196 

14.53 ± 0.68 

1560 

1.72 

0.02324 

0.583 

18.48 

2400 

0.4623 ± 0.0040 

0.0428 

7.74 ± 0.24 

1570 

1.72 

0.02324 

0.583 

18.41 

2100 

0.5256 ± 0.0049 

0.0428 

8.83 ± 0.27 

1570 

1.08 

0.00925 

0.479 

14,57 

2100 

0.4159 ± 0.0028 

0.0140 

21.4 ± 1,2 

15^0 

1.24 

0.01217 

0.510 

16.00 

2100 

0.4567 ± 0.0061 

0.0196 

16.74 ± 0.81 
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275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 
302 

303 

304 

305 

306 

307 

308 


309 

310 

312 
312 

313 
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Table 1. continued 


T 

d 

B 


m 

t 

m 


K 

K 

mm 

cm 2 


mg 

s 

mgmiri" 1 

/ 

Pa 

15’"9 

1.08 

0.00925 

0.479 

14.18 

1800 

0.4726 ± 0.0030 

0.0140 

243 ± 1.4 

15*9 

1.24 

0.01217 

0.510 

15.62 

1800 

0.5213 ± 0.0060 

0.0196 

19.16 ± 0.91 

15*9 

1.72 

0.02324 

0.583 

14.83 

1500 

0.S927 ± 0.0041 

0.0428 

9.98 ± 030 

1589 

1.08 

0.00925 

0.479 

10.72 

1200 

0.5353 ± 0.0042 

0.0140 

27.7 ± 1.6 

1589 

1.24 

0.01217 

0.510 

11.91 

1200 

0.5944 ± 0.0086 

0.0196 

21.9 ± 1.1 

1589 

1.72 

0.02324 

0.583 

16.80 

1500 

0.6729 ± 0.0077 

0.0428 

1137 ± 0.36 

1598 

1.72 

0.02324 

0.583 

18.94 

1500 

0.7555 ± 0.0091 

0.0428 

12.80 ± 0.41 

1558 

1.08 

0.00925 

0.479 

15.16 

1500 

0.6048 ± 0.0094 

0.0140 

313 ± 1.8 

1598 

1.24 

0.01217 

0.510 

16.89 

1500 

0.6743 ± 0.0080 

0.0196 

24.9 ± 1.2 

1608 

1.72 

0.02324 

0.583 

16.83 

1200 

0.837 ± 0.013 

0.0428 

14.23 ± 0.48 

1608 

1.08 

0.00925 

0.479 

13.52 

1200 

0.673 ± 0.011 

0.0140 

35.0 ± 2.1 

1608 

1.24 

0.01217 

0.510 

15.13 

1200 

0,7544 ± 0.0098 

0.0196 

28.0 ± 1.3 


“Listed in the table are the run temperature, T/K, the diameter of the cell orifice, d/ mm, the area of the effusion orifice, B/ cm 2 , the Clausing factor 
for the orifice, W& the total mass lost during the isothermal period, m/mg, the duration of the isothermal period, i/s, the computed loss rate and its 
estimated uncertainty, m/(mg-min -1 ), the factor,/, as given by eq 8 for each cell, and the value of the apparent, measured vapor pressure, P m / Pa, as 
given by eq 2 for each run. The area of the evaporating surface, A , for all cells was taken as the internal cross-sectional area of the cell (0317 cm 2 ), 
and the Clausing factor for all the cells, was 0.412. 


ward, the data from each of these cells at these specific 
temperatures was plotted in the form given by eq 9. As noted 
earlier, a linear fit to this form should give the equilibrium vapor 
pressure as the intercept and the evaporation coefficient may be 
computed from the slope. The raw data from each of the 
experimental runs are given in Table 1. These observables 
include the cell temperature, orifice diameter and mass loss rate 
given in columns 1, 2, and 7, respectively. Quantities related to 
the cell geometry are also given in the table. These include B, 
the effusion area (column 3), W& the Clausing factor for the 
effusion orifice (column 4), and the / parameter for the cell 
(column 8) as defined in eq 8 with the assumption that the 
effective evaporation area is given by the internal cross-sectional 
area of the sample cell. As described in the earlier section, the 
mass loss rate, m, and its associated uncertainty are calculated 
from the averaging procedure noted earlier. In addition, the 
length of the isothermal periods and the total mass lost during 
the time are given in columns 6 and 5, respectively, of Table 1. 
Although, the mass loss rates calculated from these starting and 
ending values should be extremely close to the values reported 
in column 7 of the table, there may be slight differences 
between the two. Several different cells were used in these 
experiments with orifice diameters ranging from (0.57 to 2.19) 
mm. The values of the measured vapor pressures, and their 
uncertainty estimates are given in column 9 of the table. 

The uncertainty in these values of P m varies based mainly on 
the uncertainty in the mass loss measurement. When 
constructing the Whitmann— Motzfeldt plot of P m versus P^J 
for the different cells, this results in a varying uncertainty in 
both the x and y coordinates of the graph. Therefore, a 
computer routine, FITEXY, was used that is capable of fitting a 
line through data when there is uncertainty in both 
coordinates. Weighting factors were constructed from the 
uncertainty, in each coordinate and given by where 

Wi = 1 /of (10) 

The code for this FITEXY subroutine was taken directly from 
ref 24. Another advantage of using this routine is that it 
provides uncertainty estimates for both the slope and 
interecept, and these were used to compute the estimated 
uncertainties in the final quantities. Plots of each isotherm were 


made, and the weighted fits were compared to unweighted fits 314 
to verify that the final data were not incorrectly influenced by 315 
the weighting factors. In approximately half of the data points, 316 
the unweighted and weighted results were essentially identical, 317 
and in the remaining fraction, the results were only modestly 318 
different. 319 

These derived values are given in Table 2. In constructing the 320 12 
Whitmann— Motzfeldt plots, the data from Table 1 reduces 321 


Table 2. Evaporation Coefficient, a f Equilibrium Vapor 
Pressure, P e q/Pa, and Enthalpy of Reaction, A fm H°(298.15 
K)/(kJ'mor^), for SiO Derived from the Measurements 
Using Three Different Knudsen Cell Orifice Sizes 


T 



A m ,H' , (298.15 K) 

K 

a 

Pa 

kJ-moT 1 

1433 

0.0089 ± 0.0040 

5.7 ± 2.2 

363.8 ± 4.7 

1443 

0.0082 ± 0.0020 

7.4 ± 1.4 

363.0 ± 23 

1453 

0.0086 ± 0.0021 

83 ± 1.6 

363.9 ± 23 

1462 

0.0086 ± 0.0021 

9.8 ± 1.9 

364.2 ± 2.4 

1472 

0.0082 ± 0.0021 

12.0 ± 2.5 

364.0 ± 2.6 

1481 

0.0085 ± 0.0025 

13.7 ± 3.2 

364.5 ± 2.9 

1491 

0.0069 ± 0.0020 

19.1 ± 4.6 

362.7 ± 3.0 

1500 

0.0070 ± 0.0019 

223 ± 5.1 

362.9 ± 2.9 

1510 

0.0068 ± 0.0019 

26.6 ± 6.1 

362.9 ± 2.9 

1520 

0.0067 ± 0.0018 

31.2 ± 7.1 

363.0 ± 2.9 

1530 

0.0066 ± 0.0018 

37.2 ± 8.4 

363.2 ± 2.9 

1540 

0.0068 ± 0.0018 

42.1 ± 9.0 

363.7 ± 2.8 

1549 

0.0064 ± 0.0017 

51 ± 12 

363.3 ± 3.0 

1560 

0.0067 ± 0.0017 

57 ± 12 

364.2 ± 2.8 

1570 

0.0063 ± 0.0028 

69 ±31 

363.8 ± 6.2 

1579 

0.0059 ± 0.0027 

82 ± 39 

363.5 ± 6.8 

1589 

0.0059 ± 0.0027 

94 ± 46 

363.9 ± 7.0 

1598 

0.00S6 ± 0.0027 

111 ± 59 

363.6 ± 7.9 

1608 

0.0054 ± 0.0027 

128 ± 74 

363.7 ± 8.8 

down to the 19 points given in 

Table 2. A minimum of three 

points were 

used for each of the final values 

given in this table. 

Shown in 

the table are the 

isothermal 

temperature, the 

evaporation coefficient, and the 

equilibrium 

vapor pressure. A 


plot of the evaporation coefficient values and their associated 326 a 


E 
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327 uncertainty is given in Figure 2 . These values range from 

328 approximately 0.005 to 0 . 009 , and there appears to be a slight 



T/K 


Figure 2. Evaporation coefficient, a, for a mixture of silicon and silica. 
Shown in the figure are experimental data from this work, O, and 
values reported by Rocabois et aL, — , and Gunther, . 

329 decrease in these values with temperature. On the basis of the 

330 relatively large uncertainty in these values, it would likely be 

331 hazardous to ascribe any significance to this trend. For this 

332 reason, it is more reasonable to report a single value of ( 0.007 

333 ± 0 . 002 ) for this temperature range. 

334 In 1958 , Gunther estimated the evaporation coefficient for 

335 the silicon/silica reaction to be approximately 41 (T 3 , and this 

336 value is shown in the graph as the dashed— dotted line over the 

337 experimental temperature range reported by Gunther. 2S The 

338 values reported in this work are of the same, low order of 

339 magnitude as those of Gunther but are slightly higher. In 1992 , 


Rocabois et al. studied both the evaporation coefficient and the 3 « 
equilibrium vapor pressure over amorphous silicon monoxide 341 
and mixtures of silicon and silica . 16 One focus of this work was 342 
to study the stability of nominally amorphous silicon monoxide. 343 
Rather than being a stable compound, research suggested that 344 
silicon monoxide exists as an intimate mixture of silicon and 345 
silica over very small domains. Data were taken by these 346 
authors using multiple Knudsen cells and mass spectrometry. 347 
Measured evaporation coefficients for the silicon/silica 348 
mixture were much smaller than those for the amorphous 349 
silicon monoxide. For vitreous silica, these values reportedly 350 
ranged from 21 (T 4 to 8 - 1 CT 3 over the temperature range (1097 351 
to 1489 ) K, while those for amorphous silicon monoxide were 352 
approximately an order of magnitude higher . 16 In addition to 353 
vitreous silica, Rocabois et al. also studied the evaporation 354 
coefficients for the reaction of silicon with cristobalite, finding 355 
even lower values of the evaporation coefficient for this 356 
reaction. In this case, they reported a constant evaporation 357 
coefficient of ( 1.1 ± 0 . 5 ) 1 CT 3 over the temperature range 358 
(1172 to 1404 ) K. Rocabois et al. provided a fit to their 359 
evaporation coefficient data for their vitreous silica and silicon 360 
data, and this fit is shown as the solid line in Figure 2 . There is 36i 
some overlap between the data range covered by Rocabois et al. 362 
and the results presented in this work, and there seems to he 363 
reasonable agreement between the magnitude of the evapo- 364 
ration coefficient in this range, although the two sets of data 365 
seem to follow very different temperature trends. 366 

A plot of the equilibrium vapor pressure values from Table 2 367 
is given in Figure 3 by the open circles. A fit to the equilibrium 368 f 3 
vapor pressure values is given by 369 

log 10 (P/Pa) = (13.25 ± 0.89) 

(17900 ± 1300) 

T/K (11) 



Figure 3. Comparison of silicon monoxide vapor pressure data over a mixture of silicon and silica. Shown in the figure are experimental data 

reported in this work, O, and a linear fit to these data, — , a fit to experimental data from Rocabois et al., , a fit to experimental data taken by 

Kubaschewski and Chart, , experimental data points from Huang et aL, 0 , and a fit to data taken by Shomikov and Archakov, •••. Also shown for 

comparison is experimental data for silicon monoxide evaporation, □, reported in a previous study and a fit to these data, ■. 


F 
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370 over the temperature range (1433 to 1608) K and is shown as 

371 the solid line in Figure 3. This fit was weighted based on similar 

372 weighting factors as given in eq 10, but only for the y-values 

373 since the estimated uncertainty in the temperatures are equal. 

374 Due to the low evaporation coefficients of the silicon— silica 

375 mixture, the uncertainties in these fit parameters are more than 

376 twice as large as the uncertainties in the fit for silicon monoxide 

377 vapor pressure over silicon monoxide powder. 7 Although there 
3^8 is considerable uncertainty in the equilibrium vapor pressure 

379 values, the linear fit to the data is very good. 

380 The equilibrium vapor pressures have also been used to 

381 calculate the third law enthalpy of reaction values, and these are 

382 given in column 4 of Table 2. In computing these values, 

383 thermodynamic data for silicon and silicon dioxide were taken 

384 from the National Institute of Standards and Technology/Joint 

385 Army, Navy, Air Force (NIST-JANAF) tables. 26 It should also 

386 be noted that these heat capacity data used to generate the free 
38 ? energy functions for the third law analysis are identical to the 

388 values used in the study of nominal SiO^) evaporation. In 

389 addition, it is assumed that the uncertainty in these enthalpy 

390 values arises entirely from the uncertainty in the vapor pressure 

391 data. These values range from (363 to 365) kjmor 1 and do 

392 not seem to follow any temperature trend. Therefore, for this 

393 work the third law enthalpy of reaction at 298.15 K, 

394 A nai H :, (298.15 K), is taken as (363.6 ±4.1) kj-mol" 1 . Using 

395 eq 11, the enthalpy of the reaction at a mean experimental 

396 temperature of 1521 K is (355 ± 25) kjmor 1 . Again, using the 

397 heat capacity data for silicon and silicon monoxide as noted 

398 above, this gives the second law enthalpy of reaction at 298.15 

399 K for the data as (355 ±25) kjmor 1 . There is considerable 

400 uncertainty in this second law value due to the large uncertainty 

401 in the slope of eq 11. These values are in fair agreement and 

402 certainly agree within the experimental uncertainty reported for 

403 each value. 

404 A comparison with some of the most recent measurements 

405 of SiO vapor pressure is also given in Figure 3. In 1974, 

406 Kubaschewski and Chart measured the vapor pressure of SiO 

407 over a mixture of silicon and silica using a thermo gravimetric 

408 and the Knudsen effusion method. 27 Data were taken in the 

409 temperature range of (1270 to 1600) K, using different orifice 

410 diameters and extracting estimates of the evaporation 

411 coefficient and equilibrium vapor pressures in the same manner 

412 used in this work. These authors fit their vapor pressure data to 

413 the following equation: 


Iog 10 (P/Pa) = 13.613 - 


17850 

T/K 


( 12 ) 


414 and this fit is shown by the dashed— dotted line in Figure 3. The 

415 slope of this fit and the current data from this work are very 

416 close, but the magnitude of the Kubashewski and Chart vapor 

417 pressures are larger than in this work. Evaporation coefficients 

418 were not specifically reported in the work of Kubaschewski and 

419 Chart but can be extracted from their vapor pressure data. 

420 Shomikov and Archakov calculated the evaporation coefficient 

421 for the silicon + cristobalite mixture as (1.43 ± 0.16) 10^ from 

422 Kubaschewski and Chart's vapor pressure measurements, 28 

423 As previously mentioned, Rocabois et aL studied both the 

424 evaporation coefficients of silicon monoxide gas over silicon 

425 monoxide powder and mixtures of silicon/silica in 1992. 

426 Although they found large differences in the evaporation 

427 coefficients for both of these systems, the vapor pressures in 

428 both cases were essentially identical within the experimental 

429 uncertainty in their work. Individual vapor pressure values were 


not reported in their work, only graphs of fits to their data. The 430 
fit to their SiO vapor pressure data is shown as the long dashed 431 
line in Figure 3. In a previous study with silicon monoxide, very 432 
close agreement was found between the measured vapor 433 
pressures and data from Rocabois et al. with an overlap in the 434 
experimental temperature range. These data for the evaporation 43 s 
of silicon monoxide powder are shown as the square data 436 
points and a fit to these points (short dashed line) in Figure 3. 437 
For the current data for silicon and silica evaporation, the vapor 438 
pressures fall somewhat below these data for silicon monoxide 439 
powder, and agreement between the two is just at the edge of 440 
the estimated uncertainty between the two measurements. 441 
Later in 2000, Shomikov and Archakov published results 442 
from a study of the evaporation of silicon monoxide, focusing 443 
on the determination of the evaporation coefficients for 444 
amorphous silicon monoxide and mixtures of silicon and 44 s 
silica. 8 Similar to Rocabois et al., these authors also used the 446 
Knudsen effusion combined with mass spectrometry. For a 447 
mixture of cristobalite and silicon, Shomikov and Archakov 44 s 
found the evaporation coefficient to be (1.65 ± 0.10)T(T 3 , a 449 
value which compares very favorably with results from 450 
Kubaschewski and Chart (1.43 ± 0.16) *1CT 3 and Rocabois et 451 
al. (1.1 ± 0.5)* 1(T 3 . 28 452 

Although the equilibrium vapor pressures for silicon 453 
monoxide were not reported by Shomikov and Archakov, 454 
they can be estimated using data presented by the authors in 455 
graphs of P m and reported / parameters for their cells. These 456 
data were digitized for the four cells used to study the 457 
evpaoration coefficients of the crystalline silica and silicon 458 
system. In constructing the Whitmann— Motzfeldt plot, the cell 459 
with the largest orifice diameter did not seem consistent with 460 
the remaining data so only the three cells with the smallest 46i 
effusion orifices were used to estimate the equilibrium vapor 462 
pressure. This estimate of the vapor pressure from Shomikov 463 
and Archakov is given as the dotted line in Figure 3. These 464 
extrapolated, equilibrium vapor pressure values are approx- 465 
imately 20 % higher than the experimentally measured vapor 466 
pressures taken with the cell with the smallest effusion orifice. 467 
These values certainly he within the experimental uncertainty of 468 
this work, yet have a slightly different slope. 469 

Huang et al. studied the vapor pressure over a silicon and 470 
silica mixture at pressures much higher than those available via 471 
Knudsen effusion. 29 The goal of this work was to estimate the 472 
equilibrium vapor pressure of SiO over this system to aid in 473 
modeling the evaporation rate of SiO in the Cz— Si system. 474 
These authors built an ampule of silica, placed an amount of 475 
silicon within the ampule, and then evacuated and sealed the 476 
ampule. The ampule was then placed under a carbon heater 477 
within a vacuum chamber charged with an atmosphere of argon 478 
gas. As the ampule was heated, the silicon and silica would exert 479 
a vapor pressure that would distort the ampule. The outside 480 
argon pressure was then adjusted to eliminate the distortion of 481 
the cell, thus matching the SiO vapor pressure within the cell. 482 
These authors measured SiO vapor pressure at three temper- 483 
atures with estimated uncertainties in these values of ± 130 Pa, 484 
and these values are shown in Figure 3 as the three points at the 48 S 
highest temperatures of the plot. Since the cell of Huang et al. 486 
was sealed, the pressure exerted on the walls of their ampule 48 : 
would equal the equilibrium vapor pressure of silicon monoxide 488 
over a mixture of silicon and silica, and the evaporation 489 
coefficients do not play a role in this system. 490 

It is interesting to note that there is reasonably good 491 
agreement on the equilibrium vapor pressure of silicon 492 
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493 monoxide between the collections of data when the evaporation 

494 coefficient is correctly applied. For example, there is reasonable 

495 agreement with the vapor pressure data over amorphous silicon 

496 monoxide powder in spite of the fact that the evaporation 

497 coefficient for this system is approximately an order of 

498 magnitude higher than for the silicon/silica system. Even in 

499 the case where there are different forms of silica used (vitreous 

500 versus crystalline), the equilibrium vapor pressures for the 

501 silicon monoxide gas are approximately equal when the 

502 evaporation coefficient is correctly taken into account. 

503 ■ CONCLUSION 

504 The evaporation coefficients and equilibrium vapor pressure of 

505 silicon monoxide gas over a mixture of silicon and vitreous silica 

506 has been studied using Knudsen effusion. The rate of mass loss 

507 of silicon monoxide gas was measured using a thermogravi- 

508 metric balance. Due to the very low evaporation coefficient for 

509 this system, several experiments with Knudsen cells of different 

5 10 effusion orifice sizes were made, and the data taken at constant 

511 temperatures were plotted and extrapolated to zero orifice size 

512 to get equilibrium vapor pressures. Vapor pressure and 

513 evaporation coefficient data were measured over the temper- 

514 ature range of (1433 to 1608) K. The resulting evaporation 
sis coefficient, (0.007 ± 0.002), was very low and in reasonable 

516 agreement with measurements by other research groups. The 

517 enthalpy of reaction at 298.15 K for this reaction was calculated 

518 via second and third law analyses as (355 ±25) kJ-moT 1 and 

519 (363.6 ± 4.1) kj-mor 1 , respectively, and was in agreement 

520 within the uncertainty in the values. It is thought that 

521 amorphous silicon monoxide is not a true comp ound, bu£ ^ 

522 rather an intimate mixture of silicon and silicor/^nonoxide?" 

523 Therefore, the kinetics of the evaporation reaction will depend 

524 upon the value of the evaporation coefficient, hut the 

525 equilibrium vapor pressure for the (nominally) amorphous 

526 silicon monoxide system and an ordinary mixture of silicon and 

527 silica should be identical. A comparison of the data for the 

528 current silicon and silica reaction with a previous experimental 

529 work with silicon monoxide powder shows that this is the case 

530 and provides further evidence that solid silicon monoxide is an 

531 intimate mixture of small domains of silicon and silica rather 

532 than a true compound 
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